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Abstract 

The Schrodinger equation is thoroughly analysed for the isotropic oscillator in the 
three-dimensional space of constant positive curvature in the spherical and cylindrical 
systems of coordinates. The expansion coefficients between the spherical and cylindrical 
bases of the oscillator are calculated. It is shown that the relevant coefficients are expressed 
through the generalised hyper geometric functions 4-F3 of the unit argument or 6j Racah 
symbols extended over their indices to the region of real values. Limiting transitions to a 
free motion and flat space are considered in detail. 

Elliptic bases of the oscillator are constructed in the form of expansion over the spherical 
and cylindrical bases. The corresponding expansion coefficients are shown to obey the 
three-term recurrence relations. 
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1 Introduction 



Starting with the classical works by Schrodinger [1] , Stivenson [3] and Infeld [2] systems with ac- 
cidental degeneracy in spaces of constant curvature have attracted attention of many researchers 
in connection with nontrivial realization of hidden symmetry in these problems and with pos- 
sible applications, especially to constructing many-particle wave functions [4], nonrelativistic 
models of quark systems [5] and solutions of the two-center problem [6]. 

Essential advances in the theory of systems with accidental degeneracy have been made by 
Nishino [7], Higgs [8], Leemon [9] and [10, 11, 12, 13]. It has been shown that the complete 
degeneracy of the spectrum of the Coulomb problem and harmonic oscillator on the three- 
dimensional sphere in the orbital and azimuthal quantum numbers is caused by an additional 
integral of motion: an analog of Runge-Lenz's vector (for the Coulomb potential) and an analog 
of Demkov's tensor (for the oscillator). However, in contrast with the flat space the integrals 
of motion for the Coulomb problem and isotropic oscillator do not form the Lie algebra as the 
relevant commutators are nonlinear. The latter does not allow one to restore, respectively, the 
algebra or the group of hidden symmetry. Later in the works [14,15] it has been shown that 
as an algebra of hidden symmetry one can use quadratic algebras of the general type, the so- 
called Racah algebras. Systems with hidden symmetry for the harmonic potential and those of 
Winternitz-Smorodinsky's type in the three-dimensional space of constant curvature were also 
studied by using the technique of path integrals in the papers by Barut, Inomata and Junker 
[16,17], Grosche [18] and Grosche et al. [19,20,21]. 

A possible way of determining a group of hidden symmetry of systems with accidental 
degeneracy is determination of the expansion coefficients between different bases obtained af- 
ter the separation of variables in the Schrodinger equation. Such interbasis expansions have 
first been considered for the "sphere-cylinder" transitions (isotropic oscillator on the sphere), 
"sphere-parabola" transitions and those between spherical and elliptical bases (for the Coulomb 
potential on the sphere and hyperboloid) in [14,15,25]. It has been shown in [19] that like for the 
Helmholtz equation [22] variables in the Schrodinger equation for the potential of the isotropic 
oscillator on the three-dimensional sphere are separated into all the six orthogonal systems of 
coordinates: spherical, cylindrical, sphero-conical, two elliptic and ellipsoidal coordinate sys- 
tems. 

The aim of the present paper is the description of solutions of the Schrodinger equation in 
the spherical, cylindrical and two elliptic systems of coordinates and the calculation of expansion 
coefficients between the corresponding bases. Note that the solution of the Schrodinger equation 
for the isotropic oscillator in the spherical system of coordinates was found in [8,9,14] and in 
the cylindrical and elliptic systems of coordinates are presented for the first time. 

The paper is organised as follows. Section 2 presents some known results related to the 
Schrodinger equation for the three-dimensional space of constant curvature. Section 3 is devoted 
to the solution of the Schrodinger equation for the potential of the isotropic oscillator in the 
spherical and cylindrical systems of coordinates. Section 4 is the calculation of coefficients of 
the interbasis expansion between spherical and cylindrical bases using the explicit expression 
for the wave functions of the isotropic oscillator. In Section 5 the elliptic bases of the isotropic 
oscillator are constructed as expansion over the spherical and cylindrical ones. 



1 



2 The Schrodinger equation and integrals of motion 



The Schrodinger equation in the space of constant curvature has the form 



(1) 



where Alb is the Laplace-Beltrami operator that in an arbitrary system of coordinates is given 
by 
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Choosing a metric of the space of constant curvature in the form (r 2 = XiXi) 
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we derive the following expression for the Laplace-Beltrami operator: 
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where R is the curvature radius 
Pi = - 
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dxi R 2 dxk 
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and the following commutation relations hold: 
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It is easily seen that within the limits of the flat space, i.e. as R — > oo, the operator Pj 
corresponds to an ordinary momentum operator, and the Laplace-Beltrami operator (2) turns 
into an ordinary Laplace operator in the flat three-dimensional space E 3 . 

The issue of generalising the problem of isotropic oscillator for the spaces of constant cur- 
vature with the use of the conformally flat metric in the classical mechanics has obviously been 
solved for the first time in H, where, in particular, an additional integral of motion character- 
istic of an oscillator interaction was found. Later, in || it has been shown that if the metric of 
the curved space is chosen in the form of (3), the role of the potential of isotropic oscillator in 
the flat space is played by 



V(r) 



Ltu 2 r 2 
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and an additional integral of motion has the form: 



1 ll 2 uj 2 
D ik = - (PiP k + P k Pi) + ^-^XiX k , (8) 
l n 



which in the limit of large R exactly transforms into Demkov's tensor [ 23] . For the operators 
Li and the following commutation relations are valid: 

[Dij, Lk] = i (eikiDji + ejkmDim) , 

[Diki Dji] = — j2 — {buLkj + SikLij + 5ijL k[ + 6jkLa) + \\Lij, D lk } 



+ {Lu, D k j} + {L k j, Du} + {L kh Dy} j , L ik — XiP k — x k Pi. 

where {, } means the anticommutator of two operators. 

The three-dimensional space of constant positive curvature can also be realised geometri- 
cally on the three-dimensional sphere S3 of the radius R, imbedded into the four- dimensional 
Euclidean space, i.e on the hypersurface 

Ql + QiQi = R 2 , 

where the coordinates change in the region < R 2 and to each value of qi correspond 
two points on the sphere. Relation between the coordinates x« in the tangent space and 
(/i = 0, 1, 2, 3) is given by 

1i = 1 — ^ == r> 9o 



l + r 2 /R 2 ^/l + r 2 /R 2 

obtained under the mapping from the center of the three-dimensional hypersphere onto the plane 
tangent to the "North pole" . Such a parametrization of the space of constant curvature is often 
called in literature the "geodesic parametrization" p4j and in a one-to-one manner reflects only 



the hemisphere (in this case the upper one) or the sphere with identified diametrically opposite 
points. 

In the coordinates g M we have 

1 nr i ( 9 d \ T . d 
P * = -R N * = R [ q *Wo ~ J ' U = ~ ieim dq~ k (9) 

and 

*™ = -± (Nf + L 2 ) 
where the operators Li Ni are generators of the group 0(4) 

[Li, Lj] = ieij k L k , [Li, Nj] = ieij k N k , [Ni, Nj] = ieij k L k 
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The potential of the isotropic oscillator is given by the symmetric function 

v ( r ) = v M = t! r 1 -q2 /R 2 > <? = <i + <i + <i (io) 

with respect to the upper and lower hemispheres, equals zero at the poles of the sphere and has 
a singularity at the equator. An additional integral of motion is given by the expression 



Dik = _l f NN + NkN .) + g*gfc (in 

2R2 K l fcT fe %J ~r ^ l _ q 2j R 2 V ) 

and like in the case of the flat space leads to separation of variables in the Schrodinger equation 
in more than one system of coordinates. 



3 Solution of the Schrodinger equation 



3.1 Spherical basis 

In the spherical system of coordinates 

qi = R sin \ sin i? cos (p, q 2 = R sin \ sin "d sin ip, 
q3 = R sin x cos qo = R cos x, 

0<X<7r, < $ < tt , < y? < 2vr, 
the oscillator potential has the form 



V = 

Choosing the wave function according to 
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where ^ m (^, is an ordinary spherical function [31]], after separation of variables in the 
Schrodinger equation we have 
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and making the substitution 
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we arrive at the equation without the first derivative of the Poschl -Teller-type 
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whose general solution is well known |pq| . The requirement of regularity of the wave function 
Z(x) at x = and tt/2 leads to quantization of the isotropic oscillator energy 
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where the principal quantum number N — 0, 1, ... is related with the radial and orbital quantum 
numbers by N = 2n r + I, and the following notation is introduced: 
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Note that the degree of degeneracy, like in the case of motion in the field of the harmonic 
isotropic oscillator in the three-dimensional Euclidean space, is equal to (N + 1)(N + 2)/2. 
The solution of the quasiradial Schrodinger equation (13), orthonormalised in the interval 

x e [o, f ] is 

z( x ) = zux) 
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where P^ a '^(x) are the Jacobi polynomials. 

Let us consider the limit of the flat space. It is easily seen that at large R [y — > XR 2 , 
X = fxto/h 2 ) the formula (14) is used to restore the formula for the energy spectrum of the 
three-dimensional oscillator 
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Transition from the spherical system of coordinates on S3 to the relevant system of coordinates 
on E 3 is accomplished within the limit R —* 00, x ~ + and x ~ r/R where r is the radius vector 



in the three-dimensional flat space | 26fl . Using the known relation for the Jacobi polynomials 



where L"(x) are the Laguerre polynomials, and taking account of the limiting relations 
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we immediately get that 
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where Rj>ii(r) is the orthonormalised spherical radial wave function of an ordinary three-dimensional 
isotropic oscillator in flat space p5| . 

The second interesting limit is the transition to a free motion. As v — > (u — > 0) we have 



]im E V N (R) = E° N (R) 



h 2 (N + l)(iV + 3) 
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Comparing the above-derived expression with the formula for the energy of a free motion of 
particles on the sphere - ^jy^ we get that J = N + 1 and J = 1, 2, and consequently, the 
ground state with J = is missing in the limiting spectrum. As for the oscillator spectrum 
(N — I) is always even, within the limit of a free motion (J — /) takes odd values and at fixed J 
there exist only states with I = J — 1, J — 3, .... and, correspondingly, the degree of degeneracy 
of the limiting spectrum is smaller than (J + l) 2 , as it should be for the free motion on the 
sphere. 



Further, using the transformation |34 

X Pln~ U) \2X 2 - 1) = 



r(A)r(m + 3/2) 

r(A + m + i)r(i/2) 



(19) 



connecting odd Gegenbauer polynomials with the Jacobi polynomials and passing from the 
quantum number iV to J, we come to the function 



2 l+1 n 



v->0 



(J+iXJ-Q ! 

v^F \| (J + l + l)\ 



(smx) l C l jt)(cosx), 



which, with an accuracy to a factor, \/2 corresponds to the solution of the free Schrodinger 
equation on the three-dimensional sphere with the impenetrable barrier at the equator (x = 
tt/2). 

3.2 Cylindrical basis 

In the cylindrical system of coordinates 



qi = Rsin a cos <fii, q 2 = Rsinasm(f) 1 , 
= R cos a sin 2 , Qq — R cos a cos 02, 



the 



< a < tt/2, < 0i < 2tt, -tt < 2 < tt, 
potential of the isotropic oscillator is written as 

...2 R 2r 
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Choosing the wave function in the form 



^(0 1 ,«,0 2 ; J R) = ^=$(a)ir(0 2 ) 



/B? V27T 

after separation of variables we arrive at two differential equations of the Poschl-Teller-type 

d 2 M / m 2 -\ A-±\ 

-—+ [£-— -a-^ ^ M=0, 21 

acr \ sin a cos^ a J 

^+ 4A -^f + ^ = ' ( 22 ) 

a0 z y sm cos z J 

where (sin a cos a)~ 1//2 M(a) = $(a) = 4f- + f , G [0, |] £ is determined by expression (14). 
The spectrum of constants is determined by: 

A = (n 3 + v + l) 2 , 5 = (2n + n 3 + \m\ + v + 2) 2 , 

where the quantum numbers n 3 and n run the values 0, 1, 2, .... Assuming the principal quantum 
number N to be equal to N = 2n+ \m\ + n 3 , we get the formula (15) for the isotropic oscillator 
energy. For the cylindrical basis we get the following expression: 

*(0!,a,02;i2) = * v Nnna {<k,<*,toR) = -7m^N\m\n 3 ^)K 3 (^)^= (23) 
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where the functions K^fo) and ^^, m , n3 (a), normalised in the interval 02 £ [— f , f], « £ [0, 
are 
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Note that to a cylindrical system of coordinates corresponds an additional integral of motion 

G?0 2 H COS 2 02 



Within large R the cylindrical system of coordinates on the sphere turns into an ordinary 
cylindrical system of coordinates (p,<f),z) in the Euclidean space E 3 ||26| . Passing to the limit 
R — > oo and a, 02 — * 0, and assuming 

sma ~ a ~ — , sin 02 ~ 02 ~ -7;, 



as well as using formula (17) and the following relation |28 



where TC n (z) are the Hermite polynomials |p8| , we get 
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Thus, formula (23) leads to the orthonormalised cylindrical basis of the isotropic oscillator in 
the flat space. 

Within the limit of a free motion v —>■ (u; — > 0) with the use of the formula [28] 
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we have 
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Assuming + 1 = |m 2 |, |m| = \mi\ and J = AT + 1, we obtain (with an accuracy to a factor 
v2) odd solutions of the Schrodinger equation for a free motion in the cylindrical system of 
coordinates. 



4 Expansion between spherical and cylindrical bases 
4.1 Calculation of the transition coefficients 

Let us write the expansion between the spherical and cylindrical bases of the isotropic oscillator 
in the form 



2V-|ro| 
n 3 =0,l 



(26) 
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where the quantum number n 3 takes even and odd values depending on parity N — \m\. 

To calculate an explicit form of the expansion coefficients W^ 3 lm (v) it is sufficient to use 
orthogonality in one of the variables for the functions entering into the cylindrical wave func- 
tion and to fix at the most appropriate point the second variable that does not participate in 
integration. Passing beforehand in the left-hand side of the expansion (26) from the spherical 
coordinates to the cylindrical ones, according to the formulae 
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Then, substituting the asymptotic formulae derived into the interbasis expansion (26), reducing 
(sina)!" 1 ! and using the orthogonality of the functions i^ 3 (0 2 ) in the interval — | < 2 < f , we 
arrive at the following integral representation for the coefficients W Nlm {v): 
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A complete solution of the problem needs calculation of the integral in formula (27). Let us 
consider separately the cases of even and odd quantum number n 3 . Separating the interval of 
integration (28) into two intervals (— f ,0) and (0, |), after the substitution in the first integral 
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02 —* —(p2 we see that the value of the integral is just doubled due to parity (I — \m\ — n 3 ). 
Then, using the well-known transformation for the Jacobi polynomials |£ 
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after the substitution x = cos 2 cj> 2 , we come to the following two table integrals for even and odd 

™3 

r(n 3 + ^ + !)(f)! 



A l{+) dA 



2 U+ 2 
i 



T(^3 + z/)(n 3 )! 



x 



;-|m|-l 



(1 - X)-^(1 + X ) v+ %P { it^ + ^(x)P& ^\x)dx, 

1 T" 



^■N\m\n 3 ( V ) 



r(n 3 + ^ + |)(^i)! 



. i— |m|+3 

2^+ 2 



^ + z/) (n 3 )! 



fl (1 - x)^(l + x) v+1 iP { ^' v+1 *\x)P^^\x)dx. 

-1 2 2 



Using the formula for integration of the two Jacobi polynomials p0| 
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p + 1, /3 + r + n + 2, r-a-n + 1 



we immediately get that -Ajvj^Ug (v) is expressed through the generalised hypergeometric func- 
tion 4 F 3 of the unit argument 
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and analogously 
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Taking account of the known symmetry property for the series 4-^3(1) of the Saalschiitz-type 

(f-bUg-b) n 



29 



4-^3 



-n, b, c, d 
e, /, 9 



-4-^3 



(f)n(g)n 

-n+b+c+d=l+e+f+g 



—n, b, e — c, e — d 
e, b — f — n + 1, b — g — n + 1 



one can easily be convinced that both the hypergeometric functions 4-^3(1) entering into A 
can be transformed to a unique form: 
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After simple transformations we finally derive the sought formula for the coefficients of the 
interbasis expansion W, 



n3 I 
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+ 1 



(29) 



Note that the expression we have derived for W^ 3 lm (u) is independent of parity of the quantum 
number n^. 



4.2 Connection with the Racah coefficients 

The interbasis expansion coefficients (29) can also be expressed through 6j, the symbols or 
Racah coefficients of the SU(2) group, extended over their indices to the region of real val- 
ues. Comparing the expression for W l Nm n ,{y) with the representation of the Racah coefficients 
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W(abed; cf) through the hypergeometric functions 4-F 3 (l) of the unit argument [[31 

A(abc)A(cde)A(aef)A(bdf) 



W(abed; cf) 



(a + b- c)\(d + e - c)!(a - / + e)\{b - f + d)\{c- a- d + f)\ 

l) ,(30) 



(a + b + d + e + l)\ f -a-b + c, -b-d + f, -a-e + f, c-d-e 

-4-^3 



(c - b - e + f)\ 6 \-a-b-d-e-l, -a + c-d + f + 1, -6 + c- e + / + l 
where A(afrc) is 

A(abc) 



(a + b- c)!(a - b + c)\(b + c - a)\ 



\ (a + 6 + c+l)! 

and taking account of the symmetry property of the Racah coefficients 

W{abed; cf) = i(-l) a+b ~ c W(abed; cf), 

e = — e — 1, d = —d — 1, f = —f — l, 
after simple calculations we get the required formula 

WnM = (-l)^ + ^v^ + V2)(n3 + v + l)W(abed; cf), (31) 

_ JV+Jmj , _ JV-Jmj-1 _ 21-1 

a ~ 4 ' 4' 4 ' 

7 _ jV-|m] 1 w 1 1 „ 7V+H 1 £ f = H3 _|_ £ 

4 ~ l ~2~ l ~4' 4 ~ r 2' J 2~ l ~2" 

Then, using the relation of orthonormalization for the Racah coefficients [[31] 



]T ^(2c + l)(2f + l)W(abed; cf)W(abed; cf) = 5 ff , 

c 

we can write the inverse expansion in the form 

TV 

l=\m\,\m\+l 

where summation over / starts with \m\ or \m\ + 1 depending on the parity of the number 
N — \m\, and the coefficients 

can be expressed through the polynomials 4-F 3 (l) with the use of the representation (30). 
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4.3 Limiting relations 

Consider limiting transitions to the flat space and free motion in the expansion coefficients 

4-3. 1 As R — > oo the generalised hypergeometric function 4^3(1) transforms into 3-^2(1) 
according to 



4-^3 



m 
2 ' 

u+1, 



n 3 -l _ .V-Z N+l 1 ,. 1 9 
2 ' 2 ' 2 "t" ^ ~r Z 

l+\m\-nz _|_ Z-| m l~ n 3 _|_ -y 



1 



f _M 
p I 2 ' 

3-^2 S ;_|_| m |_ n3 



n 3 -l 



N-l 



+ 1 



2 ' 2 
l—\m\—nz 



1 



1 



' k 2 1 ' 2 

Having made the relevant limiting transition in gamma functions, after simple algebraic trans- 
formations we get the known formula for the coefficients of the interbasis expansion between the 
spherical and cylindrical bases of the harmonic isotropic oscillator in the flat Euclidean space 



m+ \m I 



(N - \m\ - n s )\\(N + \m\ - n 3 )!! 



2'-™3 \ (JV + Z + l)!!(JV-OH(n 3 )! 



v /(2/ + l)(/+|m|)!(/-|m|)! 



N-l 



r( 



/-|m|-n 3 



p J 2 ' 2 ' 2 

Z+|m|-n 3 _|_ j\ 3 2 I l+\m\-n 3 _|_ ^ H m h n 3 



1 



2 ' V 2 1 v 2 1 ^ 2 

On the other hand, the limiting transition to the flat space can directly be traced in the formula 
(31). Indeed, using at large R the asymptotic coupling [[!]]] 

1 



W(abe + R,d + R;c,f + R) 



2R(2c+l) 



C 



c, d—e 

a, /— e; 6, d— / 



and the symmetry property of the Clebsch-Gordan coefficients [31] 



(-1) 



a+b—c Q C >1 



a,a; fe,/3' 



we obtain that 



limW/->) 



m+|m| 2i^i 2|m|-l 

'_i — 2 4 ' 4 

v x / jV+Jmj N+\m\-2n 3 _ JV-|m|-l 2n 3 -N+\m\-l 

4 ' 4 ' 4 ' 4 



4-3.2. Assuming z/ = 0, (u = 0) and passing to the quantum numbers corresponding to 
the free motion on the sphere m = mi, n 3 + 1 = |m 2 |, N = J — 1 (note that J — I is odd and 
J — |mi| — |m2| is even), we have 



lim W 



".3 1 

Nlm\ 



m 1 + \m 1 I \Tfl2 
-1) 2 -! 



\ 



( J +HH^I )j( JH™ihH»l )!(j + / + l)!!(J - /)!! 



^/(/ + l/2)(/+ |m|)!(Z- |m|)! 

p ^-|mi| — |m2| _|_ 3 ^p ^+|mi|-|rrt2| _|_ 3 s 



[mgj—l 



|m 2 |-2 



4-^3 



2 ' 2 ' 

3 1+|to|-|to 2 | _|_ 3 l-\m\ — \m 2 \ 
2' 2 2' 2 ~ t ~ 2 



J-Z-l J+Z+3 
2 ' 2 

3 



■1)2+ 2 



v /(/ + l/2)|m 2 |^(a6erf; c/), (32) 



where 
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J-\m\-2 _ 21-1 

4 ' C 4 ' 

J+\m\ — 1 n |m2| — 1 

4 ' J ~ 2 

Let us mention an interesting fact that within the limit of a free motion in formula (32) for the 
interbasis coefficients W^] m (v) instead of the hypergeometric function 3 F 2 of the unit argument 
we have the function 4 F 3 ; and instead of the Clebsch-Gordan coefficients of the SU(2) group, the 
Racah coefficients for one fourth values of the SU(1,1) group momentum. Analogous formulae 
arose in calculating the coefficients of transition between different hyperspherical systems of 
coordinates (in the formalism of "trees") and have been analysed in ||32|| . In our case, this fact 
allows an alternative calculation of the integral A l N , m , n3 (i/) in formula (28) at v — 0. 

5 Elliptic bases 

The oblate elliptic system of coordinates ( known as the elliptic-cylindrical I) has the form 



a 
d 



j+H-i 

4 ' 

J—\m\ 



b 

e 



qi = i?sn(/i, k)dn(v, k') cos 0, q 2 = Rsn(fi, k)dn(i/, k') sin 0, 
q 3 = i?cn(/i, k)cn(v, k'), qo = Rdn(fi, k)sn(i/, k'). 

-K<n<K, -2K' <u< 2K', < (p < 2vr, 

where the elliptic Jacobi functions of the variables a and f3 have the moduli k and k', respectively, 
k 2 + k' 2 = 1, and K and K' are the complete elliptic integrals. 

For the potential V in the elliptic system of coordinates we have 



V(fi, u) = -Muj 2 R 2 
Choosing the wave function \1/ in the form 



dn /zsnV 



(33) 



V,(p) = ^i(/i)^ 2 (^) 



m e Z, 



(34) 



after the separation of variables in the Schrodinger equation (|33"D we arrive at two ordinary 
differential equations 



d 2 ipi cn/idn/i dip\ 



dfi 2 



d 2 jj 2 
dv 2 



+ 



sn/i dfi 



( 2MER 2 M 2 u 2 R v 
+ 



k 2 sn 2 fx + 



V ^ 2 



Tr 



m 

sn 2 /i 



M 2 u 2 R A k' 2 



Tr 



dn 2 /i 



,12 



snz/cni/ dip2 
dnz/ dv 



+ 



(2MER 2 M 2 u 2 R 4 \ , 2 

+ r9 dnV + 



-X q (k; R)ipi, 

k 2 m 2 M 2 u 2 R 4 1 



(35) 



V n 2 



Tr 



dn 2 v 



h 2 sn 2 v 



+X q (k;R)ij 2 



(36) 
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where the quantum number q enumerates the elliptic separation constant X q (k;R). Excluding 
from the equations fl35f ) and (|36|) energy E, we come to the following operator 



A 



&; 2 sn 2 /i — dn v 



2 d 2 2 2 92 . cn/idn/i 2 5 , a , 2 snz/cnz/ 9 
dn ^7—5- + rsn /i— H dn 1/— k k 



dn 2 di* 2 sn/i dn dnz/ <9z/ 



+ 



M 2 cu 2 R 4 dn 2 u + k 2 sn 2 fi - 1 dn 2 z/ + Psn 2 /i <9 2 



7i dn /isnV dn z/sn 2 /i cfy? 2 

;i _ A; 2 )L 2 - fc 2 /^ + k 2 h\ + (k' 2 )^^ + k 2 ^^H, (37) 

n n 



whose eigenvalues are X q (k; -R), and eigenf unctions are given by expression (|34]). Let us introduce 
a new operator according to 

1 f» „ l2 ,M 2 uj 2 R 4 , 2 2MR 2 TT , 2r2 l 

= L 2 -aR 2 D 33 , (38) 

where a = & 2 /(l — fc 2 ) G [0, 00). As is known |2"7| , a transition from the oblate to the prolate 
elliptic system of coordinates on the sphere can be obtained under the transformation k — > 
ik/k', k' — > 1/k'. In this case, the operator fl38|) for the prolate system of coordinates has the 
form 

Q^ 1 - = L 2 + k 2 R 2 D 33 . (39) 

The latter formula allows one to describe both the elliptic systems of coordinates uniquely with 
the use of the operator 

3 = L 2 - aR 2 D 33 , (40) 

where a G [— l,oo). For positive a we have the oblate system; and for a G [—1,0], the prolate 
elliptic system of coordinates. 



5.1 Expansion of elliptic bases over the hyperspherical and cylin- 
drical ones 

Thus, we have seen in the previous sections that all three oscillator bases: hyperspherical, 
cylindrical and both the elliptic ones, are eigenfunctions of three complete sets of operators 
{H, L 2 } L 2 Z }, {H, D 33 , L 2 } and {H, 3, L 2 } so that 

L 2 q Nlm ( X ,$,v) = l(l + i)* mm {x,#,'p), (41) 
D 33 ^ Nn . Am (ipi,a,Lp 2 ) = (n 3 + v + 1) 2 ^/ Nn3m (ip 1 ,a,ip 2 ), (42) 
Q^ n9m (/i, p, ip) = X q (a;R)^ nqm (n,v,Lp). (43) 

The operator equations (]4"l|)-(f43"l) allow us to construct elliptic bases of the isotropic oscillator 
on the sphere as a superposition over the hyperspherical and cylindrical bases. 
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Now, let us write the sought expansions: 

N 



* N9m (ji,v,(p) = T l Nqm (a-R)^ Nlm {x,$,y), (44) 

Z=|m|,|m|+1 
AT— |m| 

^NqmiPi = UNqmi^R^Nmmi^Pua,^). (45) 

n 3 =0,l 



Consider the expansion (f4~4"|). Substituting fl4"3| ) into the operator equation (^2[), we find 

1 N 
__ + i) _ Ag (a; J2)} T^ m = £ ^(AbV, (46) 

where 

(AhV = / ** Nlm (D 33 )y Nl , m( m. (47) 

To calculate the integral (flTf ) we use the expansion of the spherical basis over the cylindrical 
one and equation (42) for the eigenfunctions of the operator D 33 . As a result, we come to the 
following expression for (D 33 )w: 

N—\m\ 



{D w )„,= £ W n N \ m W n N ], m {n 3 + v+l)\ (48) 
Then, using the three-term recurrence relations for the Racah coefficients pi[ 



^'r^^M^ri+^W^H (49) 



B c = J(a + b + c + 2)(-a + b + c+l)(a-b + c+l)(a + b-c) (50) 



where 



x yj(d - I + c + l)(d + I - c)(d + I + c + 2)(-d + I + c+ 1), 
A c = [a(a + l)-b(b + l)][d{d+l)-l(l + l)] + c(c + l)[a{a + l) (51) 
+ b(b + 1) + d(d + 1) + /(/ + 1) - c(c + 1)] - 2c(c + + 1) 

and the orthogonality property 

iV-]m| 

E W^ m W^ m = 6 u ,, (52) 

"3=0,1 



we have 



(Ab)«' - - (2Z _ 1)(2Z + 1) W + " (2/ + 1)(2/ + 3) W 

16 



where 



B l = —J(l - \m\ + - \m\ +2){l+ \m\ + 1)(Z + \m\ + 2){N + I + 3){N - I) 
16 v 



x yJ(N + l + 2u + 4)(N -l + 2u + l), (53) 

d = i(4(A^+l)(A^ + 3) + 2(2|m| 2 -l)+4z/(2A^ + 2z/ + 5)-(2/-l)(2/ + 3) 
8 I 

(4H 2 -l)(2iV + 3)(2iV + 5 + z/) ) 

(2/-l)(2/ + 3) J 1 ] 

Substituting the matrix element ( |53"1) into (^), we finally arrive at the three-term recurrence 
relation 

(2/-l)(2/+l) ^ 2T ^ 2m(a; R) + {^ [Kl + Xq{a; R)] ~ C ) T ^ m(a; R) 

16 

+ (2/ + l)(2 i + 3) B ' r ^ (a;fl ^°- < 55) 
^Vgm( a i -R) = T N g m (a; R) = 0. 

for the expansion coefficients T l Nqm (a; R). The recurrence relation fl55l) is a system of homoge- 
neous equations which has to be solved with the normalization condition 

N 

£ \T l Ngm (a;R)\ 2 = l. 

i=|m|,|m|+l 

Eigenvalues of the elliptic separation constant X q (a; R) are calculated from the condition for the 
determinant of the system of homogeneous equations to be equal to zero (j55|). 

Consider now the expansion of the elliptic basis over the cylindrical one. In a similar 
way, like in the calculation of the expansion coefficients fl44]), we get 



JV-|m| 

{\{k-R) + aR 2 {n, + v + lf}u n N \ m {k-R)= £ U# qm (k; R)(L 2 ) n3n ,, (56) 



n' 3 =0,l 



where 

(L 2 )n 3 n> 3 = J y* Nn3m L 2 y Nn , 3m( m. (57) 

The integral in fl5"T|) can be calculated if one uses the expansion of the cylindrical basis over the 
spherical one and then the symmetry property of the Racah coefficients 



(5* 



a,b,c 1 _ (a,l,f 
d,l,fj Id, 6, c 

and the three-term recurrence relation (|49|). As a result of simple calculations we obtain the 
expression {U^\ m = U n3 ) 

B n3 U n3+2 + [C n3 - X q (k; R) - aR 2 {n 3 + v + l) 2 } U n3 + B n3 _ 2 U n . 3 . 2 = 0, (59) 
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where 



1 



(n 3 + 2u + 2)(n 3 + 2)(n 3 + l)(n 3 + 2i/ + 3){N + \m\ + n 3 + 2v + 4) 



(n 3 + 1/ + l)(n 3 + z/ + 2) 2 (n 3 + i/ + 3) 
x \J(N + \m\ - n 3 )(N - \m\ - n 3 )(N - \m\ + n 3 + 2v + 4), 



U(iV + 2) 2 



2 I 



+ u(2N + 2i/ + 5) + (H 2 - 2) - (n 3 + u)(n 3 + v + 2) 



1/(1/ + 1)(AT + \m\ + v + 2) (AT - \m\ + z/) 
(n 3 + z/)(n 3 + z/ + 2) 

As in the previous case, a homogeneous system of equations should be solved together with the 
normalization condition 

N-\m\ 
"3=0,1 

and the separation constant can again be determined from the condition of equality to zero of 
the corresponding determinant of the homogeneous system of equations (|59|). 

In conclusion, we would like to mention that in the limit of the free motion v — > 0, the 
three-term recurrence relations (55) and (59) transform into those for the expansion coefficients 
of the elliptic basis over cylindrical and spherical ones, which have been obtained in the work 
[27]. 



6 Conclusion 

In the present paper, a first step has been made to completely study the Schrodinger equation 
and to calculate the coefficients of various interbasis expansions for the potential of the isotropic 
oscillator on the three-dimensional sphere in different orthogonal systems of coordinates. We 
have calculated interbasis expansions for the "sphere-cylinder" transition and also constructed 
solutions of the Schrodinger equation in both the elliptic systems of coordinates as expansion 
over the spherical and cylindrical bases. In contrast with the "sphere-cylinder" transition, 
in which the transformation coefficients are expressed through the generalised hypergeometric 
functions 4F3 of the unit argument or through the Racah coefficients, extended over their indices 
into the region of arbitrary real values, the transition coefficients for the elliptic bases are defined 
by the three-term recurrence relations and cannot be written down explicitly. 

In this paper we have considered the sphero-conic and ellipsoidal bases of the isotropic 
oscillator. In separating variables in the Schrodinger equation for the sphero-conic system of 
coordinates we arrive at the quasiradial equation (13) considered in Sect. 3 and two standard 
Lame equations that are derived in separating variables for the Helmholtz equation on the 
two-dimensional sphere. As a result, the solution of the Schrodinger equation is 
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where the function Z v Nl [x) is determined by expression (16), and the explicit form of the Lame 
polynomials Au(a) can be found in |28], |36 |. 



For the ellipsoidal basis of the isotropic oscillator, after separation of variables in the 
Schrodinger equation we obtain three identical equations containing two ellipsoidal separa- 
tion constants Ai, A2. This means that in contrast with the cases considered we deal with the 
two-parametric spectral problem. Application of the method of constructing solutions of the 
Schrodinger equation as expansion over simpler bases leads to two many-termed recurrence rela- 
tions determined by the cubic matrix. Obviously, the simplest way of constructing an ellipsoidal 
basis consists in applying the Niven method |37j allowing one to write down a solution to the 



Schrodinger equation in terms of zeroes of the wave function and reduce the problem to a system 
of relevant nonlinear equations. This kind of investigation is beyond the scope of the present 
paper and will be carried out elsewhere. 

In conclusion, we would like to thank L.G. Mardoyan and V.M. Ter-Antonyan for useful 
interesting discussions. 
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